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Abstract 

We study the interaction of a fast moving particle in the Quark Gluon Plasma with linearized 
hydrodynamics. We derive the linearized hydrodynamic equations on top of an expanding fireball, 
and detail the solutions for a static medium. There are two modes far from the jet ~ a sound 
mode and a diffusion mode. The diffusion mode is localized in a narrow wake behind the jet 
while the sound mode propagates at the Mach angle, cos{9m) = c^/c. A general argument shows 
that the strength of the diffusion mode relative to the sound mode is directly proportional to the 
entropy produced by the jet-medium interaction. This argument does not rely on the linearized 
approximation and the assumption of local thermal equilibrium close to the jet. With this insight 
we calculate the spectrum of secondaries associated with the fast moving particle. If the energy 
loss is large and the jet-medium interaction does not produce significant entropy, the flow at the 
Mach angle can be observed in the associated spectrum. However, the shape of associated spectra 
is quite fragile and sensitive to many of the inputs of the calculation. 
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I. INTRODUCTION 



One of the major findings at RHIC is jet quenching, the suppression of high transverse 
momentum particles jfj . This strong suppression is attributed to the energy loss of partons 
traversing the dense medium formed in a high energy heavy ion collision. Different micro- 
scopic mechanisms for the jet energy loss have been proposed 0, 0, 0, IE 0] (sse Ref. for 
a review). 

Ultimately the energy "lost" by these fast partons is shared among the constituents of 
the medium. In a previous paper we argued that the fate of the deposited energy and 
momentum should be described by hydrodynamics^ 0. Hydrodynamic models describe 
the RHIC data reasonably well and motivated this suggestion. The suggestion was also 
motivated by the experimental observation of particles associated with a high px trigger 
at a particular azimuthal angle, A0 = vr — 1.2 rad, relative to the trigger particle [l^]l5l|. 
Within linearized hydrodynamics, jets induce flow fields with a conical shape similar to 
the supersonic flow past planes. It was argued that this conical flow could provide an 
explanation for the observed correlation. Discussion of alternative explanations to the large 
Pt correlations at RHIC can be found in Refs. SEi 11113. 



In this first paper we solved linearized hydrodynamics in a static homogeneous medium, 
and found that two modes can be excited - a sound mode and diffusion mode 0. The 
relative strength of these modes depends on the jet-medium interaction and will be clarified 
in this work. Subsequently, the effect of the fireball expansion was studied |2fl|], and an 



attempt to incorporate more realistic mediums has been made [21|. Further, the nonlinear 
hydrodjTiamic response to the jet was estimated with 2+1 dimensional hydro code Jl^]. The 
conical flow was not observed in the azimuthal correlations of that work; possible reasons 
for this are discussed in section IVIII Finally, a parton transport calculation based on the 
AMPT model also found large angle di-hadron correlations [22^, .2^1 . 

In a second paper two of us studied the conical flow in an expanding fireball with a 
variable speed of sound and found that the amplitude of the sound wave increases jis^ . 
Furthermore, we have shown that a first order phase transition with a vanishing speed of 
sound, would lead to more peaks in the di-hadron correlation function. These additional 
peaks have not been observed. 

In this paper we present a systematic study of the conical flow within the linearized 
theory. In section IH] we start by providing an overview and explaining the approximations. 
Then in section UTTl we review our calculation for a static medium and find the two modes of 
linearized hydrodynamics jsj. In section HVi we describe potential flow, first for an ideal fluid, 
and then with viscosity. Section HVl also provides the linearized equations which propagate 
small perturbations over a general hydrodynamical expansion. In section we relate the 
energy and momentum loss of the jet to the outgoing sound waves at large distances. This 
result ties the initial energy loss to the normalization of the final spectrum of secondaries 
that is calculated in section IVIllI Our work in section |3 also delineates two microscopic 
models for the jet-medium interaction based on entropy production. These two models are 
described in sect ions IVIl and IVTTl and their corresponding correlation functions are calculated 
in section IVIIII Finally, many technical derivations are delegated into appendices. 



see Ref. for a similar idea by Stocker. 
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FIG. 1: A schematic picture of the flow created by a jet moving through the firebah. The trigger 
jet is moving to the right away from the origination point (the black circle at point B). Sound 
waves start propagating as spherical waves (the dashed circle) from the origination point. The 
companion quenched jet is moving to the left creating a wake of matter (shaded area) and adding 
to the sound wave. The head of the jet is a non-equilibrium gluonic shower formed by the original 
hard parton (black dot A). The solid arrow indicates the flow velocity which is perpendicular to 
the shock cone at the angle 6m, cos{9m) = Cg/c ~ 0.55. 

II. OVERVIEW 

A schematic picture of di-jet production in a nucleus-nucleus collision is given in Fig. [H 
For simplicity we will only consider a di-jet pair with both particles at mid-rapidity in a 
central Au-Au collision. One jet serves as the trigger particle and is biased toward the 
edge of the nucleus due to the strong jet quenching. The other jet then travels a distance 
2-Rau ~ 10 fm and is either completely or partly absorbed. Experimentally a companion jet 
with pt ~ 10 GeV seems to be absorbed while a jet with larger transverse momentum seems 
to "punch through" and re-appear at an angle A0 = n relative to the trigger jet. 

In general the rapidity of the two jets are not identical and this will complicate the 
experimental interpretation of our results. The particular configuration we are considering 
can be selected with three particle correlations and currently there is an experimental effort 
in this direction. 

Far from the jet we can calculate the correlation of T'^'^ with the passage of the jet using 
linearized hydrodynamics. Generally there are two modes - a sound mode and a diffusion 
mode. The diffusion mode is concentrated in a narrow wake behind the jet while the sound 
mode propagates forward at the Mach angle, cos(6'a/) = c^/c ~ 0.55. The total momentum 
loss can be related to the amplitude of the wake and the amplitude of the sound wave. Under 
the reasonable assumption that the total energy and momentum loss of the jet are equal we 
reach the stronger conclusion which relates the entropy production of the jet to the relative 
strength of the two modes. 

The equations for the sound wave are given by Eqs. ()4.2|) . and ()5.3p . The unknown 
function in this formula, dF/ dx, is related to the momentum transfered to the sound wave - 
see Eq. (16. 5p . Similarly the flow flelds in the wake are given by Eq. (j5.4j) , and the amplitude 
A in this formula is flxed by the momentum transfered to the wake - see Eq. (j5.13|) . This 
momentum transfer is in turn related by Eq. (j5.14j) to the total entropy produced by the 
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jet. In summary, by specifying the total rate of energy loss and entropy production the flow 
fields at large distances are determined. 

We next estimate the medium modifications at the head of the jet where the jet loses 
energy through collisional and radiative processes. Since the jet energy is large compared 
to the typical scale of the medium, the jet acts as a point source. The energy deposited by 
this source is distributed into the surrounding fluid through highly dissipative processes. To 
estimate the order of magnitude of the initial modification, consider a jet that losses certain 
amount of energy per unit length, dE/dx. This energy is absorbed by the medium over a 
dissipative length scale of the order of the sound attenuation length, Tg = Arj/Sw, where rj is 
the shear viscosity and w = e + p is enthalpy of the medium. Thus we compare the energy 
deposited by the jet in this typical length with the energy of the fluid in a deposition 
volume of r^. 



To estimate this initial modification, we use the formula from Ref. for radiative energy 
loss 

-^ = ^^LlnA, (2.2) 

dX 8 Xg Xg ^ ^ 

2 

where ^ = g is the transport coefficient. Perturbative estimates for this parameter give 

q = O.e'GeVVfm for T = 200 MeV 0. Setting Xg ^ 1 L = 5fm, and Cr = 3 

(gluons), we obtain, dE/dx ^2.7 GeV/fm. Similarly, perturbative estimates for the sound 
attenuation length give with = 1/2, Tg ~ 0.18/T [ll|, [l^- The conjectured lower bound 
for Ts is I/SttT [28|. For the energy density of the unperturbed medium (e) we will take the 
QGP value as measured on the lattice, e ~ 12T'^ jl^]. For T = 200 MeV, we finally conclude 
that 

/? ^ X F 

1^ ^ '^^ " ^ 36 - 100 > 1 . (2.3) 

^fluid e X F3 

(The range of values is set by the range in F^.) Thus, this quantity is numerically large 
although it is suppressed suppressed by in perturbation theory. 

Since this number is greater than one, the jet is surrounded by its own small "fireball" (of 
size Fs or more) where variation of the thermodynamic quantities is very large and hydrody- 
namics cannot be applicable. Outside of this region, there is a domain where gradients are 
small enough that viscous hydrodynamics can in principle be used, but the behavior of the 
fluid is non-linear, dissipative, and possibly turbulent. We will not discuss these complex 
regions in the present work. 

Our current objective is to study what happens far from the jet, where the situation 
becomes less violent and more tractable. Specifically, it is the region where the flow velocity 
and pressure modifications may be considered small compared the unperturbed medium, 
i.e. the region where linearized hydrodynamics can be used. We note that it is possible to 
calculate the drag force on an airplane using linearized hydrodynamics, by computing the 
momentum flow through a large cylinder around the plane. The analogous computation for 
jets is performed in section |Vl 

In summary one may separate space-time into three regions 



4 



1. A region very close to the jet, i.e. a region of size Tg, where strong dissipative effects 
happen. 

2. A region close to the jet where hydrodynamics is still nonlinear and possibly turbulent. 

3. A region far from the jet where the perturbations and gradients are small. In this 
region one may use linearized hydrodynamics and the motion is laminar. 

Since the primordial interaction of the jet with the fluid is complicated and non-linear, 
we will only study flows in the third linearized region, and impose only general constraints 
on the total energy, momentum, and entropy deposited by the jet. 



III. LINEARIZED HYDRODYNAMICS IN A STATIC MEDIUM 

In this section we will review our results of Ref. |3| where we studied small perturbations of 
a homogeneous relativistic baryon free fluid at rest. The linearized hydrodynamic equations 
are written in terms of the small quantities 

e = (5r°° , g' = 6T^' . (3.1) 

The remaining components can be expressed in terms of these quantities by means of the 
equation of state. For a fluid with shear viscosity rj and vanishing bulk viscosity these 
remaining components are 

w ^ ' 

where = dp/de is the speed of sound of the fluid, w the enthalpy and {d'^g^) is symmetric 
traceless tensor 

(9V> = 9V + 5V-^^W. (3.3) 

By writing the energy and momentum conservation equations dfj,5T^^ = 0, the linearized 
hydro equations can be written after a spatial Fourier transforms as 

dte + ikgi, = , 

dtgi + tcik e + ^^1— k^gL = , (3.4) 
3 eo + 

dtgT + -^k^gT = , (3.5) 
Co + Po 

where g = (^l^ + St- Thus the equations of linear hydro decouple into two different modes 
- a sound mode (Eq. ()3.4|1 ) and a diffuson mode (Eq. ()3.5|) ). The excitation of these two 
modes depends on the initial conditions. 

In order to study these initial condition, we considered the disturbance produced by an 
infinitesimal displacement of a particle moving along the x direction with velocity v in a 
time interval dto at time to- According to the axial symmetry of the problem, the most 
general expression for the initial disturbance is 

(^dto{t = to,x) = eo(x,r) (3.6) 
gdto {t = to,^)= 9o{x, r)(5" + Vgi{x, r) 
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The source functions eo{x,r) and gi{x,r) excite only the sound mode, while the remaining 
function go{x,r) excites the diffuson mode. The particular value of these functions depends 
on the interaction of the jet and the fluid in the near region. As argued in the introduction, 
it is difficult to find the exact functional form for these sources and they are only constrained 
by the total energy, momentum, and entropy deposited by the jet. 

The hydrodynamic disturbance due to these infinitesimal displacements are calculated as 
an integral of the sources 

rO'^(x) = j {K>: (x - y) eo(y) + K>^^ (x - y) go{y) + K^^^ (x - y) gM) , (3-7) 

where the three kernels K^, Kj^^ and Kj^^ are 0. 

= E (i^*^'.(r.).a'^'^'.(r.)) . (3.8) 



i=± 

Here we have defined the functions 



+ I 0, d'd.^-£p^ + 6'^ 3 ,^ I , (3.9) 

^ — u \ * / 



1 (r±c,t)2 3 1 Z±- 

P±{rs) = -^^e-^^ , P(-r.) = e ^ . (3.11) 

In what follows we will try to clarify further the meaning of the different modes as well as 
their excitations. We will start addressing the problem in the case of inviscid fiuid and we 
will show how to connect to the viscous case. The results presented for dihadron azimuthal 
distribution will be performed through these solutions. 



IV. RELATIVISTIC POTENTIAL FLOW 

In nonrelativistic laminar fiow it is customary to introduce a potential, so that the velocity 
can be expressed as a gradient, v = V0. Such fiows are clearly irrotational as the V x v 
vanishes. The conservation of circulation for inviscid fiuids guaranties that potential fiow 
remains irrotational throughout the evolution of the liquid. Sound waves, in particular 
correspond to small potential disturbances over some given fiow. 

One can define a relativistic analogue of the potential fiow for an ideal fiuid in which the 
vorticity of the fiuid vanishes . Let us first recall the definition of the the vorticity in the 
relativistic case jl^l 

uJ^r = di^u^^KK, (4.1) 
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where h^'^ = rj'^'^ — -u^m^ is the projector into the space perpendicular to the velocity. From 
this definition, it is simple to see that any fiow of the form 

fu, = d,^, (4.2) 

with / some function, will be irrotational. The scalar field 0, called the potential, is the 
relativistic analogue of the nonrelativistic potential fiow. 

Thus, one can try to find solutions of the ideal relativistic hydrodynamics of the form 
Eq. ()4.2j) . In fact, the function / can be determined by requiring that the projection of the 
energy momentum conservation equation into the space perpendicular to u'^ vanishes, 

Kpd^T'-' = . (4.3) 

In the case of an ideal baryon free fiuid, these equations can be expressed as, 

s {u^d, (Tu,,) - d,,T) = 0, (4.4) 

where s and T are the entropy density and temperature of the fiuid. It is clear then that, 
by setting / = T a velocity field of the form of Eq. ()4.2|) is a solution of these equations 
j29|. Thus, imposing irrotationality of the fiow leaves us only one equation to solve for the 
baryon free fiuid, the entropy continuity equation. 

9^(sM^) = . (4.5) 

This equation can be used to find a non-linear equation for the potential (see appendix 
lA which is equivalent to the ideal equations of motion of an irrotational fiuid. 

Let us conclude by some examples of potential fiow. Obviously, the static homogeneous 
baryon free fiuid (which we will use in the rest of the paper) is potential. A slightly less 
trivial example is the boost invariant Bjorken solution j3l|. In appendix lA 21 we show how 
to find this solution from the potential. Finally, as in the nonrelativistic case, sound waves 
are irrotational disturbances of the hydrodynamic solution and are therefore potential. 



A. Small perturbations of a potential flow 

We now will study small perturbations over a given potential hydrodynamic solution. 
Even though the interaction of the jet with the background fiuid leads to very violent fiows 
in the near zone, the small perturbation approximation is always valid far enough away from 
the high energy particle. 

In the region where the perturbation is small, we will have modifications both from the 
thermodynamic quantities and from the velocity field, T', u'^. 

Let us start with the case where the perturbation of the fiuid does not change the ir- 
rotational character. In this case, the modified hydrodynamic fields can be described by a 
perturbation of the potential ^ 

+ (4.6) 



^ Here and in the rest of the paper, the modification of the aU the hydrodynamic fields but the potential 
will be denoted by '. The modification of the potential denoted by tp. 
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Expanding the definition of tlie potential Eq. ()4.2|) to first order we find 

T' = u^d^^ , (4.7) 

with h^i, = Tj'^^ — u^u^ . Note that since {u + u'Y = 1, we have m • m' = to first order in the 
perturbation. 

In a general case, since the interaction with the high energy particle is very complicated, 
the modified field may not be irrotational. However, one can still define a modified potential 
ip that describes part of the perturbation. In order to do so, we can separate the perturbed 
velocity field u'^ = Ufj, + Rfj,, such that and T' are related to the potential ip as in Eq. ()4.7p . 
The rotational component R is then the part of the perturbation that cannot be expressed 
as Eq. (j4.7|) and, thus, is linearly independent from f/; as f/ is orthogonal to u, so is R. 

We now discuss the equations of motion for the perturbations. As we have defined U 
such that it can be thought of as coming from a potential, the perturbed transverse equation 
Eq. (|4.3|) depends only on R (and the background quantities) 

T R'd^d.cl) + d''<pd,iT R^) = 0, (4.8) 

with (j) the potential of the background medium. The remaining equation is the the entropy 
continuity equation for the perturbation 

d^^su^")' = . (4.9) 

Using the definitions of the speed of sound for a baryon free fluid = (s/T) dT/ds and the 
perturbations Eq. ()4.7|) we can re-express Eq. ()4.9|) as a second order differential equation 
for if: 

Let us mention that such second order differential equation can be expressed as the 
d'Alambertian of a scalar field in a particular metric jS^ that is dependent on the ther- 
modynamic and velocity fields of the background medium: 

-^d, (v^G^-c},0) = -^d, is R^ , (4.11) 
G'"' = Cs- fr^^^ + (4 - iKw' 



where the rotational component acts as sound source. The reduction of the sound equation 
to a Klein-Gordon equation in a gravitational background is well known in nonrelativistic 
fluids where G'^'^ is known as the acoustic metric 0^. 

Let us finally write these general equation for a simple case, the static homogeneous 
medium. In this case, the equations of motion for the rotational and irrotational field 
decouple. For the irrotational field, we find the standard wave equation 

\d^<P-V'<P = 0. (4.12) 



^ This emission (that vanishes for the static case) may be important in the final observation of the conical 
flow 
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For the rotational mode, we obtain a non propagating equation 

doRi = 0. (4.13) 

These two modes coincide with what we called sound and diffuson mode in in the case 
where the viscosity of the medium t] vanishes. As also found there, only the sound mode 
propagates, while the energy /momentum of the diffuson mode remains localized close to the 
position of their deposition. Thus, as far away from the source of the perturbation (the jet) 
the modification of the hydrodynamical fields are small, we expect that only potential flow 
should be found in this region, since perturbation can only arrive to long distances form 
the source due to propagation. Thus the irrotational flow should be concentrated in the 
region close to the jet path. This is also what happens in the case of flow past bodies in 
nonrelativistic fluids. 



B. Viscosity and Potential Flow 

Our previous discussion of the potential flow was made for inviscid fluids only, while now 
we will study the effect of the shear viscosity. We show below that in the case of static and 
homogeneous background medium, one can still find a potential solution that describes the 
system. For non-relativistic fluids this fact has been shown in js^]- 

The energy momentum tensor for a liquid with non- vanishing shear viscosity is |29l | 

T^"" = wu^u'' -r]^"'p + r]{W) , (4.14) 
where w is the enthalpy, p is the pressure, = h^^d^ , and 

3 

In the case of a viscous fluid, the entropy equation is obtained by projecting u^d^T^'^ and 
leads to 

d,{su^') = ^{V^u^)d,u^. (4.16) 

For a static or homogeneous medium, there is no entropy production due to viscosity. The 
defining feature of the potential flow is zero vorticity, i.e. the existence of solutions of the 
hydrodynamic equation in the form of Eq. (j4.2j) . For small perturbations on top a static fluid, 
viscosity will modify the function / in this equation. As in the inviscid case, the functional 
form of / is determined by requiring that the perpendicular hydro equations Eq. ()4.3|1 hold 
for the viscous case 

In order to find / for the small perturbation, let us expand 

f = T + f = T + T' + m (4.17) 

The first order terms in this expansion coincide with inviscid fluid since viscous corrections 
vanish for the static homogeneous medium. If there were no viscous corrections, /' would 
coincide with the perturbed temperature T' as in the ideal case. Thus the function m 
depends on the viscosity. 
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Now consider a small perturbation of the hydrodynamic fields, T' and u'^. As in the 
previous section, we can split the velocity perturbation in a potential part that does not 
have vorticity If^ and a rotational part that carries the vorticity, W^. We define now a 
potential related to the irrotational part of the velocity perturbation using Eq. (j4.17j) as 

T' + m = uf^d^if , (4.18) 

This definition ensures that if is zero, the perturbed field does not introduce vorticity. 

After this identification, the perturbation of the perpendicular equation Eq. ()4.3p includ- 
ing shear viscosity can be written in terms of and m only as 

doRi + d,m = hsV^R, + r,5,VV , (4.19) 

where the right hand side comes from the viscous tensor which has two independent contri- 
butions coming from the irrotational and the rotational part of the perturbed velocity field. 
Thus, identifying 

m = r,VV , (4.20) 

we absorb all the viscous effects coming from potential flow into a modification of the 
function /. After finding this modification due to the viscosity, we can use the definition of 
the potential Eq. (j4.18j) into the (perturbed) entropy continuity equation to find an equation 
for if 

- c^VV - r,5tVV = . (4.21) 
The equation for the irrotational part is 

doR^ = ^TsV^Ri. (4.22) 

Eq. 1)4.2111 is the wave equation including dissipation coming from viscosity. This equation 
for the potential coincides with the equation for the energy density that one would find from 
our analysis in [8], what is expected as the potential and the energy density in this case 
are related by a time derivative. Eq. ()4.22|) is the diffusion equation for the rotational fluid 
which coincides also with the equation for the diffuson in Q. 



V. ENERGY AND MOMENTUM OF SMALL PERTURBATIONS 

As argued in the introduction, an accurate matching of the initial disturbance to the 
subsequent flow field is not possible. We can only constrain global quantities such as the total 
energy and momentum deposited. In this section we study the total energy and momentum 
of the modified fields in order to learn about the possible mechanisms of interaction. 

To this end, we consider a very high energy particle that travels with finite velocity v 
though an infinite homogeneous ideal baryon free fluid. We also assume that the motion 
has taken place forever and is therefore static in the rest frame of the particle. The analysis 
here parallels the discussion of non-relativistic supersonic flow past finite bodies [29^ . 
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As shown in section IIV A[ the equation for small potential perturbations of the potential 
flow can be readily computed from the acoustic metric, leading to the usual wave equation 
We will assume that the potential depends on the combination x + vt, with v the velocity 
of the jet. With this constraint, the wave equation far from the disturbance reads 

P'dl<P-dl^ = 0. (5.1) 

Here = x + vt, ± denotes the coordinates transverse to the jet trajectory, and 

In the region where this equation applies, we can write general solution in the following 
form: 



27rs v^(x + vt- e)2 - 

The function dF/dx{^) characterizes the source and will be elucidated in section IVIl 

From the discussion in section HV Al we know that the rotational part of the flow field 
(in the small perturbation regime) does not propagate and its evolution is determined by 
the viscosity. Examining the kernel in Eq. 13.91 the rotational field far from the fluid is 

^ (2vr|;r.(x + vt)) {'{lT,{x + vt)) ] ' ^^'^^ 

where the constant A will be determined later^. This field also depends on x + vt. 

In order to proceed further, it is beneficial to go to the rest frame of the jet, where 
the unperturbed medium moves toward the jet with velocity v. In the jet rest frame, the 
fields do not depend on the proper time coordinate r and depend only on the longitudinal 
coordinate ^. x = 1 + '^t) ■ 

The flow picture is sketched in figure |21 The non hydro core co-moving with the jet has 
a typical transverse size a. The rotational flow remains concentrated along x to the right of 
the jet in a narrow wake with a transverse extent that grows as ^^/TsX/^'y. By contrast, the 
sonic disturbance propagates and at sufficiently large distances from the jet is concentrated 
at the Mach angle within a typical size given by a. Note that for both types of excitations 
the modified fields are always small sufficiently far away from the jet. The fiuid to the left 
of the jet cannot be modified by the high energy particle and fiows with a steady velocity, 

V. 

In this stationary situation, the energy and momentum loss can be calculated by com- 
puting the momentum and energy fiux out of a large cylindrical region surrounding the 
jet. Disregarding the left end cap at x = ~oo (where there is no modification), the surface 
integral includes the cylinder surface C, and the right end cap, S. The irrotational fiuid 



^ The same expression is found in non relativistic fluids past bodies where such field describes the non- 
turbulent wake behind the body [i^ 

^ Note that this conclusion depends on the fact that the velocity of the jet does not change significantly on 
the passage thought the medium which should be a good approximation for a high energy particle. 
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FIG. 2: Sketch of the flow picture in the fluid jet rest frame. The non-hydro core (solid region) 
serves as a source for the hydrodynamic fields. The irrotational or potential component of the 
velocity field (wavy region) propagates out of the source and leads to a small disturbance at large 
distances. The irrotational part(circled region) remains concentrated along the jet axis within a 
transverse size that grows as ^Jx- 



only affects a small subregion Sq of the right end cap, S of area ~ Y sxl^l- The potential 
flow contributes to cylinder C and to the rest of the end cap, E — Sq. Thus the energy lost 
by the jet is given by 



dr 



I {T'P)'+ I (T0^)'+ I [T^^y (5.5) 



where ' denotes the variation with respect of the unperturbed field. As in the region So the 
main modification of the fields is due to the rotational disturbances, the two first terms in 
Eq. ()5.5|) are the energy loss due to the potential flow, while the last is due to the rotational 
flow, which we denote as —dE^/dr. Since far away from the jet the linear approximation is 
valid, we have 

(r°^)' = [sT^\^y = T7 {s^v'^y + S7t;^ (T7)' . (5.6) 

In both these expression we can identify the first term as the modified entropy flux times 
a constant factor (T7). The second term vanishes in the potential region since the fluid is 
stationary and we have (^7)' = d-rf^ = 0. Next note the entropy produced per unit time is 

"^^1- f (s^^P)' + f (s^^x)' + ^ (5.7) 



dr Jc is-So dr 

where dSf/dr is the entropy carried by the rotational flow in the Sq region. The two 
quantities dSf/dr and —dEf/dr are not independent (what we show at the end of the 
section) but we prefer to keep then explicitly in what is next to make our argument clearer. 
Putting these results together the energy loss is 

dE. _ dS^ dEf ^ dSf ^ ^ 
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Thus, in the rest frame of the jet, the energy loss of the jet is related to the entropy 
production minus corrections due to the rotational flow. 

Since the high energy particle remains almost on shell as it transverses the medium and 
looses energy, we conclude that the energy and momentum loss rate are linked by the particle 
velocity v as dE/dt = vdP/dt. In the particle rest frame we conclude that 

dEi 

—^^0. (5.9) 
dr 

This relation immediately allows us to arrive to an important conclusion about the interac- 
tion mechanism. From Eq. ()5.8|) and Eq. ()5.9p we conclude that if the interaction mechanism 
produces significant entropy, then rotational flow is needed and the importance of the ro- 
tational flow with respect to the potential flow is related to the contribution of entropy 
production to the energy loss. Thus, we can distinguish two different types of interactions: 

1. Isentropic interactions. The interactions of the jet and the fluid are such that no 
entropy is produced in the process. In this case the energy /momentum deposited can 
be calculated from the the far field sonic wave and it is quadratic in the perturbation. 

2. Non isentropic interactions. The main mechanism of energy deposition proceeds by 
transferring of heat into the fluid, creating new entropy. As a consequence, rotational 
flow is needed, which is concentrated in a narrow wake behind the fluid. 

Before continuing we perform a similar analysis for the momentum loss. Integrating over 
a large cylinder as before, the momentum loss is 

f (T^py+ I (T«)'+ / (r«)' . (5.10) 

df Jc JE-So Jt.0 

Now, in the region of small perturbation we find 

(T'^py = (v'^TP^)' = V {TP^y + v'^TP^ (5.11) 

where we identified a term proportional to the energy loss (due to the potential part of the 
flow). Denoting, d [P^) /dr the momentum loss due to the rotational flow (which contains 
the sum of the momentum flux through the surface Eq and vdEj^/dr), the total momentum 
loss is given by 

dr dr Jc is-E„ dr 

where we separated once again the contributions for the rotational and irrotational part of 
the flow as in Eq. ()5.8|) . This equation will be useful in the next section. 

To complete this section, let us relate the energy and momentum loss and the entropy 
rate in the wake to the amplitude of the rotatinal fleld. Boosting Eq. ()5.4|) to the jet rest 
frame we flnd = and is a Gaussian in the transverse coordinate p with width 

depending on x- From Eq. (j5.1H) 

d(P^)^ r 
^ ^ ' -= -Tsv-f^ / = -sTvA (5.13) 



dr 



So 
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In the same way, by expanding Eq. ()5.6|) and using the fact that dEj/dr = we obtain 



7 (ir 7 (ir T^y^ dr 

VI. ISENTROPIC INTERACTIONS 

As argued in the previous section, a purely potential excitation of the fluid by a high 
energy particle implies that no significant entropy can be produced in the interaction. In 
this case, the perturbed flow field in the far region is given by Eq. ()5.3|) and depends on a 
unknown function dF/ dx{^) that characterizes the source. To clarify the physical meaning of 
this source, let dF/dx{0) mark the smooth beginning of the source which moves with the jet 
velocity. Further assume that the typical scale of variation of this function A is much larger 
that the transverse distance at which Eq. (|5.3p is applicable, i.e. the source is elongated. In 
this case, for p ^ A and far away from the tip of the source, we can approximate Eq. (|5.3|) 
with 

-dF/dx{x + vt)\n( . . ) , (6.1) 



27rs ' ' ' \2{x + vt) 

Thus, for transverse distances short compared to the length of the jet, the velocity of the 
perturbation f ^ decreases as 1/p, leading to finite radial flux of entropy out of a cylinder 
close to the jet. With this normalization dF/dx{^) is the entropy flux per unit length at a 
position (^) with respect to the beginning of the jet. This is completely analogous with the 
non-relativistic case, where the coefficient is related with the mass flux ® |29l] . 

As already stated, this analysis is identical to the non-relativistic case of flow past a solid 
body. As in that case, we see that there is the formation of sound waves, or vorticity free 
excitation that fall slowly far form the body. This slow fall means that the waves carry 
energy and momentum out of the body, and we can compute the momentum drag from the 
far field solution where the linearized approximation is justified. This loss by sonic emission 
is called sonic drag. 

In order to compute the momentum losses and to use standard techniques from text books 
29l |. it is beneficial to go to the rest frame of the jet by performing a Lorentz transformation. 



as done in the previous section. We will see that in the final expression for the energy and 
momentum loss all the 7 factors drop out, so we can use the expressions obtained for the 
case of finite mass for massless particles. 

As is a scalar, transforming to the jet rest frame means simply the transformation 
of the coordinates. Thus, by performing the change of variables ^ — > ^7 and defining 
A{^) = dF/dx{C,/'y), we obtain 

^ -M, (6.2) 



^ In the nonrelativistic case, matter is pushed out of the cyhnder by the profile of the sohd body that moves 
in the medium. In our case there is no solid wall which prevents the liquid from reaching center of the 
jet. However, fluid is pushed out this region by interaction of the medium with the high px particle. We 
will not speculate further about these microscopic details. 



14 



with X the longitudinal coordinate in the jet frame and P = As we note in the previous 
section, there is no dependence of the proper time r (stationary situation). 

We now express the energy and momentum loss expressions Eq. ()5.8p and Eq. ()5.12|) for 
the particular case where there is neither entropy production nor rotational flow. In this 
case the energy loss in the jet rest frame vanishes identically. By using the definition of the 
perturbed potential, Tu'^ = d^ip and noting that the perturbed velocity field vanishes as 
X ^ oo we conclude that 



dr 

'a. 



, (6.3) 



where denotes the integral around a cylindrical surface far from the fluid. Eq. (jfj.Hj] gives 
the correct expression of the energy momentum loss by the high pt particle if there is neither 
entropy nor vorticity production. 

This is exactly the same expression that is obtained in the nonrelativistic case, but 
replacing the corresponding quantities as already explained. We can then just refer to 
standard textbooks [29] to write the final value of the integral. 



dP^ 1 T '•^ 

dr Att s 



j ^dirdi2Mii)Mi2) (In 16 - 61 - ln(4L)) , (6.4) 



where L is the longitudinal size of a long cylinder {L ^ A) and A denotes the derivative. 
The last term in Eq. (|6.4p diverges in the limit L — oo unless the flux of entropy out of the 
cylinder vanish far away from the jet. 

We now perform a Lorentz transformation to the fluid frame. We start by noting that 
according to our definition of A{S,) = dF/dx{^ /'-/), we can change the variable of integration 
in Eq. ()6.4|) and in the limit L —>■ oo the expression is independent of 7. The transformation 
dt/dr = 7 also cancels the 7 factor of the Lorentz transformation leading to an expression 
for the energy and momentum loss that is applicable when f — > c 

dE __ dP'^ _ _ dP^ 

dt dt dr ' 

dP? 1 T f°° d d 

' diidi2-rdF/dx{ii)—dF/dx{i2)\n\ii-i2\. (6.5) 



dr Atx s J _^ d^ d^ 

Note that, as a consequence of the first equation in Eq. ()6.5|1 . the wave drag is such that 
the particle remains on shell. The relative sign between the energy and momentum loss is 
due to the fact that, according to our conventions, the jet is moving in the —x direction. As 
the particle remains on shell, the isentropic assumption is a fully consistent way to describe 
the interaction of the jet and the medium. Note also that we assume that the speed of 
the propagation did not change, what should not be important for particles moving with 
velocity close to c. 



VII. PRODUCTION OF ENTROPY 

We discuss now the case where the significant entropy is produced in the interaction of the 
jet with the medium. In this case the disturbance remains close to the interaction region in 
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a narrow wake. In the wake, hydrodynamics and hnearized hydrodynamics are not justified 
until the viscosity dissipates the strength of this disturbance. At the end of the section we 
will also briefiy discuss the spectrum of particles produced for such excitations. 

In the language of infinitesimal disturbances of section II I II this case can be accounted 
for by initial condition such that eo(to)X) has a finite integral over space. In the linearized 
region the total field is found by summing these infinitesimal disturbances. In this setting, it 
is evident that the energy loss is fist order in the perturbation as, to this order, it coincides 
with the integral of the function eo 

As seen in the section |Vl the requirement that the particle remains almost on shell on 
its propagation implies that the production of entropy should be accompanied by vortex 
fiow. In the language of section UTTl it corresponds to a function go(to;x) which cannot be 
expressed as a gradient (V x v 7^ 0) and with finite integral of the gg component. Such a 
case was studied in with the following particular expressions for the initial conditions: 



€o{to,x) = go{to,x) = — — -3 , (7.1) 



where rj(to) = (^o, —to, 0, 0) is the jet position. In this case, the energy and momentum loss 
is dE/dx, and it is first order in the perturbation. Let us note that, as after setting the value 
of the fields at to we only solve the field equations for t > to, the final solution obtained as 
the addition of all the disturbances at different times to takes the form 



' dtoe{t-toW-{t,) , (7.2) 

where 5T^'^ is a solution of the linearized field equations with initial conditions specified at 
to fsection llII|) . Thus, the differential equation that our solution satisfies is, in fact, 

d^T^"" = , (7.3) 

where J depends on the fields at each to which in the particular case used is 

J' = ^^^:r^T^ ■ (1' -1' 0' 0) • (7-4) 

dx {2'n:a^) 



Let us remark here that the source Eq. ()7.3|) used in jsf is the same used in j22|] to 
solve hydrodynamic equations in the nonlinearregime in an expanding, boost invariant, 
ideal, baryon free fiuid ^. Thus, our solution is the linearized version of this problem for a 
homogeneous medium. 

After reexpresing the fields in this way, we want to connect the description of the fiuid in 
terms of section HVl and show that excitations where entropy is produced the modification of 
the fields are mainly concentrated in the near region. This is clearly seen from the rotational 
component (excited by g^^). As demonstrated in section HVAI for a static medium the 
vorticity of the fiuid does not propagate, and the disturbance remains close to the deposition 
region, i. e. in a region of typical size a the source size, and thus far from the source size the 



''' the actual calculation in [i^ integrates several of these sources along the space rapidity 77 in order to keep 
boost invariance but the conclusions should not change 
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fluid is potential. This was also found in [q] by solving the equations of motions including 
viscosity. 

Thus, if the interaction jet liquid is such that significant total entropy is produced, then 
the modification of the hydro fields are mainly produced in the near region (the 'non- 
hydrodynamic core) and, as argued in the introduction, they should not be described in 
terms of ideal hydrodynamics; one instead would need to use dissipative hydro in order to 
gain some understanding from this region. 

To conclude this section let us briefly discuss the particle production for medium exci- 
tations that produced entropy, which will be extensively done in the next section for the 
entropy free case. As found in the excitation of the diffusion mode (vorticity) leads to 
matter moving preferentially in the direction of the jet. This flow is due to the region close 
(distances of the order of a) to the jet path and leads to a complete shadowing of any Mach 
cone signal in the dihadron azimuthal distributions, resulting in correlations only at = vr. 
Similar conclusions have been reached in 22| , confirming our result for the linearized static 



theory in the non linear expanding regime. Thus, this kind of excitations, if present, leads 
to a fill up of the Mach cone in the two particle correlations associated to the jet. For more 
details we refer the reader to our short paper and to i^]. Let us remark once again that 
our calculated spectrum in Q is dominated by the near region, where the non hydrodynamic 
core is, and thus, the description in terms of hydrodynamics is not justified. 



VIII. SPECTRA AND CORRELATIONS ASSOCIATED WITH THE JET 

We now calculate the spectrum of particles associated with a jet through the induced 
disturbance of the medium. We will consider, as in the previous discussion, a jet propagating 
in a static homogeneous baryon free fluid. The jet is created at some time to = and is 
absorbed due to the energy loss at some later time tj. We take the life time of the jet tj as a 
parameter, which means that, for a fixed value of the energy loss, we are varying the energy 
of the jet. In a realistic simulation, this time would be given by the energy loss and the 
spectrum of produced jets. Finally, the system freezes out at a time tj where we calculate 
the spectrum though the Cooper Fry prescription [s^. We will always consider tj < tf so 
that the final spectrum is not dominated by our arbitrary source. 

We will only study the case where there is no entropy production, and consequently the 
energy loss is second order. We will consider the superposition of infinitesimal disturbances 
described in section IIIII and demand that no entropy is produced. This can be achieved by 
considering infinitesimal displacements such that the modification of the fields at a time to 
in a time interval dtQ is 

edto = , (8.1) 

The source size a is the characteristic distance from the jet such that the linearized approxi- 
mation is valid. The amplitude of the disturbance is C/a^ with C a dimensionless constant. 
The factor l/cx^ is introduced ad hoc such that the disturbances have the correct dimensions 
keeping C dimensionless. 

As explained in the introduction, the process by which the energy of the jet is thermalized 
is highly dissipative and we cannot describe it. Thus, we cannot fix from first principles the 
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two parameters of the initial perturbation C and a unless a good matching mechanism is 
found. Thus we will use the energy and momentum loss to relate their values. Following 
appendix El we can relate this disturbance with a far field solution in the form Eq. ()5.3p . 
After appendix ini we identify the function dF/dx in term of the parameters C and a as 

where we have introduced for simplicity a = va/cg. With this identification we can evaluate 
the energy loss using Eq. leading to 

^-^^^ (8 3) 

As expected, the energy loss is quadratic in the amplitude of the disturbance, C. The 
strong dependence on the typical size a is partly artificial and comes from our choice of 
using a as the dimensionfull parameter while keeping C dimensionless. 

We will also consider that the static medium has a non vanishing shear viscosity that 
leads to a reduction of the amplitude of the sound waves as it propagates out of the distur- 
bance. As this attenuation of the wave can be understood [i^ as a result of the dissipation 
of the mechanical energy (given by Eq. (|8.3|) ). we will use the previous formula for the 
energy /momentum loss also in the viscous case. 

Let us summarize here the parameters that we have in our calculation. The source has 
two parameters, the amplitude C and the source size a. This two can be constrained by 
the value of the energy loss. The timings are characterized by the jet lifetime tj and the 
observation time tf. In a realistic simulation, the lifetime is determined from the jet energy 
and energy loss; the freeze out, time from the hydrodynamical simulations at RHIC. Finally, 
the medium is characterized by the speed of sound c^, the sound attenuation length and 
the temperature of the static medium. We will fix the speed of sound to be = l/VS and 
we will not consider its change here in spite of our findings in |25| that showed that those 
changes are very important for the final observation and the position of the peaks. The value 
of Tg is taken close to its minimal bound F^ = I/AttT [2^. Finally, everything is calculated 
in units of the temperature, that should be of the order of the critical temperature. From 
now on we will assume that the jet moves at the speed of light v = 1. Let us finally give a 
numerical value by substituting the QGP equation of state as measured by lattice calculation 
113 (e 12T4) 

0.63^^2 ij. (8.4) 

Once we have listed all the parameters that enter in our calculation we can calculate the 
spectrum of particles induced by the jet. As the medium is static, we use as the Cooper Fry 
prescription for equal time freeze out surface as in 

dN r d'V _E^s ,0 

e ^ (8.5) 



d'^p Jy 2tt^ 

where V is the volume of the fireball, T is the temperature at freeze out and S the pertur- 
bation of the fluid. For 5 = we obtain the unperturbed spectrum from the static medium. 
When a jet goes through the fluid, 6 is related to the modified fields as: 

E 6T pv , , 

6 = -— + Y- (8-6) 
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From these expression we expect non trivial features of the pt dependence of the corre- 
lation function. As discussed in fsj, we can distinguish two different regimes in the particle 
production: 

1. Low energy particles ~ T. In this region we can expand the exponential in Eq. ()8.5|l 
and express the spectrum in terms of the energy and momentum deposited 

(Pp (27r)3 y ^ T w T w j ' ^ ' 

The term proportional to the fireball volume V corresponds to the uncorrelated con- 
tribution that is subtracted from the dihadron correlation function. The other two 
terms are due to the disturbances. As seen in Eq. (|8.7|) . the soft particles are insen- 
sitive to the particular shape of the flow field and the angular dependence is just a 
cosine of the relative angle between the observed particle and the jet. Thus, linearized 
hydrodynamics is unable to produce peaks for very soft particles. 

2. High energy particles E ^ T. In this region the large parameter E/T compensates the 
small flow velocity. The integral is then dominated by the maximum of the exponent. 
Thus, the final spectrum reflects the shape of the sonic disturbance and only points 
of maximum modification of the hydrodynamic fields contribute to the integral. This 
fact is also responsible for the absence of the conical flow when significant entropy is 
produced. In this case the region of maximum modification is close to the jet and 
the spectrum is dominated by the near field region where the applicability of the 
hydrodynamic approach is uncertain. 

Finally, we calculate the spectrum of correlated particles by subtracting the unperturbed 
spectrum (5 = 0) from the jet induced spectrum (5 7^ 0) 



dN dN 



dydcj) dydcj) 



dN 



dydd) 



(8.8) 



A. Comparison with experimental data 

From the previous discussion of the spectrum we expect that for sufficiently large energy 
particles, the azimuthal distribution of particles associated with the jet should reflect the 
modification of the fiow fields due to the particle passage. As argued in pi] this translates into 
the appearance of large angle correlations in the azimuthal dihadron distribution associated 
with the jet. Such large angle correlations have been observed by both the PHENIX ^ 
and STAR 3 collaborations. 



In this section we will fix a set of values for our parameters that give similar dihadron 
correlation than those in with no other justification as to reproduce qualitative features 
of observed on experiment. Those values are a = 0.75/T, = O.l/T, dE/dx = G.ST"^, 
tj = 8/T and tf = 10/T. The corresponding correlations functions are shown in the lower 
left panel of figure El Let us note that for a typical temperature of the plasma of the order 
of T 200 MeV, the required value of the energy loss is very large dE/dx = 12.6 GeV/fm. 
If we however would use a more realistic value of the energy loss dE/dx = 2 GeV/fm we 
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FIG. 3: Left: Associate yield dependence on associate pT for fixed source size a = 0.75/T , 
viscosity = O.l/T, tj = 8/T, tf = 10/T, and energy loss, dE/dx = lOT^ (top) and dE/dx = 
63T^ (bottom). The label values for dE/dx correspond to T = 200 MeV. The three curves are 
for IT < pt < 5T (solid), 5T < pt < lOT (dotted), (3x) lOT < pt < 15T (dashed), (lOx) 
15T < pt < 20T (dashed-dotted). (in the upper panel all the curves are rescaled further up 
by a factor 10). No large angle correlation is observed for dE/dx = WT^. For dE/dx = 63T^ 
the position of the peak shifts toward vr for lower px- Right: Experimental dihadron azimuthal 
distributions from STAR (top) Q and PHENIX (bottom) Q 

are not able to generate any large angle correlation or magnitude of the correlation for the 
intervals of momentum considered, as shown in the upper left panel of figure El * 

Even though the necessary value of the energy loss within our static approximation is 
very large, we will try to argue in the next section, where we discuss the dependence of 
the calculated spectrum on the different parameters, that expansion effects may help to 
reduce the necessary energy loss. That is why we want to point out here some qualitative 
agreement of our calculation (with large energy loss) with the experimental data. This is 
the dependence of the correlation. 

In the lower left panel of figuTe^{dE / dx = 12.6 GeV/fm) we show the pt dependence of 
the azimuthal dihadron correlation for different Pt windows. The four curves correspond to 
Pt intervals (in units of temperature) of 1-5,5-10,10-15,15-20 for sohd, dotted, dashed and 



Let us note however that, according to [2a| RHIC data seem to support values of q that are up to 4 times 
larger than those of the peturbative estimates from ^26j , which lead to energy losses that are much closer 
to our required value. 
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dashed dotted respectively. Note that we have rescaled up the two highest pt bins by factors 
3 and 10 respectively. This figure shows clearly the behavior of the correlation function. As 
claimed, at low pt, when the ratio pr/T is small (solid line), the formation of conical fiow 
in the Mach direction is not refiected in particle spectra; we obtain, in fact, that particles 
are produced mainly in the direction opposite to the jet. As px increases, the correlation 
develops, moving the peak to the left toward the Mach angle and making the width of the 
peak at this angle smaller. In the upper left panel, for smaller values of dE/dx (where all 
the curves are further rescale up by 10), the correlation always peaks at A0 = tt for the 
intervals of pt considered. This is because the fields in this case are too week. 

In figure El (right panel) we show the dihadron correlation function as presented by STAR 
(top) and PHENIX (boton) jl^]. Even though our calculation cannot be directly com- 
pared to the experimental situation, as we have not included the effect of the expansion 
(20I |2]| nor the variable speed of sound ji^ , that may be very important for the final obser- 
vation of the effect, some qualitative features seem to be reproduced. In fact, the correlation 
function shown by STAR, seems to be a broad peak with a small peak at a finite angle 
away form vr (marked by the red arrows). This could be explained from our discussion of 
the Pt dependence as in the experimental conditions of STAR, where the associated pt is 
dominated by soft particles (0.15 < px < 4GeV/c), and thus, according to the previous 
discussion, the strength of the correlation should be small, as is due to the high-pt tail of 
the associated particles. On the contrary, in the case of PHENIX, where the associated 
particles are harder (1 < p^- < 2.5GeV/c), the correlation function shows a well defined 
peak off TT. What is more, for the particular values of the parameters chosen, the amplitude 
of our calculated correlation function corresponds roughly with the experimental magnitude 
observed. 

Let us remark here that the px dependence of the correlations coming from the conical 
flow, which seems to be in agreement with the experimental situation, is c omp letely different 
from those radiative based mechanism to explain large angle correlation |l7l isi 19| . These 



model predict that the correlations should shift toward A(/) = vr as the px of the associated 
particle increases, which is the opposite to what we find. 

The position of the peak in our correlation function is at ^ tt — arccos (l/v^) = 
2.2 rad. and is set by the value of the speed of sound in the medium, which in our case 
is the ideal QGP value of = 1/3. This is of course not the case in the experimental 
condition at RHIC, where as the medium expands and cools, the speed of sound changes 
from the previous QGP value to ~ in the mixed phase and c\ ~ 0.2 in the hadron gas. 
From hydrodynamical simulations of Au-Au collisions at RHIC leads to similar (proper) 
time extension of the three previously mentioned phases of approximately r ~ 4 — 5/m. 
As argued in the Mach angle depends only on the ratio of distance traveled by the 
disturbances to the one traveled by the jet; thus the position of the peak is modified due to 
the expansion and can be estimated from the average speed of sound along the evolution of 
the medium 

cos^M = — I Csdt ^ 0.333 ^ A0 ^ TT - 1.2 rad. , (8.9) 
cr Jo 

This value agrees nicely with the experimental value shown in the experimental correlation 
function as seen in figure El b) 

Finally, in figure |31 we show the spectrum of correlated particles in the away side jet. We 
defined the away side as particles as Acf) > 1 as in We compare it with the spectrum 
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FIG. 4: Spectrum of associated particles in the away side Ac/) > 1 for a = 0.75/T, dE/dx = 142T^, 
Ts = 0.1/T (solid). Spectrum of uncorrelated particles (rescaled down by a factor 100). The 
associated yield is much harder than the inclusive due to the boosted liquid induced by the jet. 

obtained for uncorrelated particles, as if no jet were produced. In order to be able to show 
both in the same plot, we rescale down this last spectrum (that we will call inclusive) by a 
factor hundred. We observe that the correlated spectrum is harder than the inclusive one. 
This is due to the fact that the liquid is boosted due to the non zero velocity fields in the 
Mach direction. The experimental spectrum shown in fl5^ also seems slightly harder than 
the inclusive one, however the effect is by no means so prominent. On the other hand our 
inclusive spectrum is much steeper than those in jl^] as our background fluid lacks radial 
flow. Thus, even though there is no quantitative agreement between our calculation and 
STAR data in this issue, it could be due to the static approximation that we took for the 
liquid, that is definitely not true. 



B. Dependence on the parameters 

In the previous subsection we saw that in order to produce similar correlation functions 
as the experimental ones we need too large values for the energy loss. In this section we 
want to show that our calculation is very sensitive to the parameters chosen. 

In figure El we study the dependence on the energy loss (or the source amplitude) for the 
rest of the parameters fixed. We vary the energy loss by a factor of two from the central 
value of dE/dx = GST^ (z.e. 25.2, 12.6, 6.3 GeV/fm for T = 200 MeV). We observe an 
extremely strong dependence on the amplitude (note that the solid curved is rescaled down 
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FIG. 5: Associate yield dependence on energy loss for fixed source size a = 0.75/T, Ts = 0.1/T, 
tj = 8/T, tf = lO/T and lOT < pt < 20T. The three curves are for {l/20x)dE/dx = 126T^, 
dE/dx = 63T^, dE/dx = 35. 5T^ for solid, dotted and dashed respectively. 

by a factor 20). This is easily understood as from Eq. ()8.5|) we see that the amplitude enters 
in the exponent and thus small changes lead to large variations on the final observation. This 
strong dependence is the basis for our claim that expansion effects may be very important 
for the final observation of the effect on experimental conditions as it can reduce the needed 
values of dE/dx. 

In ji^ two of us found, by studying a simple model for a dynamic medium (that of a 
static fluid in an expanding universe), that the expansion of the liquid and dropping of the 
speed of sound leads to an enhancement of the ratio v/T. Estimates made for RHIC tell us 
that for disturbances originated at early times, such enhancement can be as big as a factor 
3. What this means for the spectrum is that at the time of freeze out, the amplitude of 
the disturbance can be up to a factor 3 smaller than in the static case to reproduce the 
same correlation (as the smaller amplitude is compensated by the dropping density and 
speed of sound of the background medium). Such a smaller value of the amplitude enters 
quadratically into the energy density what means that the energy loss can be up to a factor 
9 smaller than those showed in the static case. 

In figure El we study the source size dependence of the number of associated particles to 
the jet dN/dydA(j) in an interval of transverse momentum lOTy < px < 20Tj for the mid 
rapidity region y=0 for fixed values of the energy loss {dE/dx = 63T^) and the jet energy 
(tj = 8/T) at a freeze out time oft/ = 10/T. We observe that the smaller the source size the 
smaller the signal in the azimuthal dihadron correlations. This dependence is due mainly to 
the viscosity. When the source size is comparable to the sound attenuation length, viscous 



23 



^ 0.45 



■a 




0.5 1 1.5 2 2.5 3 

A(t)(rad) 



FIG. 6: Associate yield dependence on source size a, for fixed energy loss dE/dx = 63r2, = 
O.l/T, tj = 8/T, tf = lO/T and lOT < pt < 20T. The three curves are for a = l/T, a = 0.75/r 
and a = 0.5/T for solid, dashed and dotted respectively. 

effects reduce the amplitude very fast. For larger sources, the reduction of the amplitude is 
delayed till times of the order of 

tr ^ a^/Ts . (8.10) 

We also notice that as the size grows, the peak in the correlation shifts slightly toward vr. 
This is because the interference of sound waves leads to the Mach angle for propagation and 
observation times larger than the typical size of the object (as an extreme example consider 
a large object that moves an infinitesimal distance at supersonic velocity; it is clear then 
that it only can generate spherical waves). 

In figure [7| we show the dependence on the sound attenuation length F^ for cr = 0.75/T 
and dE/dx = 63T^ for three different values of the viscosity T Tg = 0.1,0.2,0.5 for solid, 
dotted and dashed hues respectively (note that the two last ones are rescaled up by a 
factor 10 and 20 respectively). As in the case of the amplitude, the correlation function is 
very sensitive to the value of the viscosity. By changing the viscosity by a factor of 2 the 
whole correlation structure almost disappears and gets reduced by an order of magnitude. 
Increasing the viscosity by a factor 5 leads to hardly no signal in the dihadron azimuthal 
distribution. This is because the dissipative effects induce a reduction of the amplitude of 
the fields. 

In figure IHl we study the dependence on the energy of the jet though the jet life time for 
a fixed energy loss. We also fixed the source size a and the viscosity F^. We again observe 
a strong dependence of the correlation strength on the energy of the jet. This is in a sense 
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FIG. 7: Associate yield dependence on viscosity for fixed source size a = 0.75/T and energy loss 
, dE/dx = 63T2, tj = 8/T, tf = lO/T and lOT < pt < 20T. The three curves are for = O.l/T, 
(5x) = 0.2/r, (20 x) = 0.5/r for solid, dotted and dashed respectively 

due to the dependence on the viscosity. As the jet has less energy it gets absorbed earher 
in the medium and thus the sonic disturbances have to propagate longer before freeze out. 
Thus, the attenuation of the waves leads to very strong suppression of the correlation signal. 
The strong dependence on the path length complicates the interpretation of the large angle 
correlations due to the conical flow, as it seems to filter out a tight range of jet energies, 
reducing the total amplitude. However this dependence on the path may not be directly 
extrapolate to the real expanding case, as in such situation, the energy loss and the viscosity 
will depend on time. Thus, even though longer times are required to observe a particle that 
is absorbed earlier, the large values of the amplitude and the smaller value of the viscosity 
at early times leads to an enhancement on the final observation. Thus, detailed studies in a 
real case should be performed before conclusions on the path length dependence. 

By the previous studies of the different parameters we want to show that the appearance 
of conical flow seems to be very fragile, given the large changes observed in the correlation 
for small changes of the different characteristic scales lead to completely different correlation 
strengths and shapes. However, if the interpretation of the large angle correlation as coming 
from conical flow survives further tests, this strong dependence reveal that this prove seems 
to be very sensitive to the details of the dissipative processes in the medium. This strong 
dependence leads us to have some hope to use this prove in order to learn about the transport 
properties on the medium. The situation is however unfortunate at the moment due to the 
freedom we have in the parameters. If we could constrain the source size by some microscopic 
model, the energy loss would constrain the amplitude of the wave, leaving us with a very 
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FIG. 8: Associate yield dependence on the jet life time tj for fixed source size a = 0.75/r energy 
loss , dE/dx = 63r2, viscosity = O.l/T, tj = 8/T, tj = lO/T and lOT < pt < 20T. The three 
curves are for tj = 8/T, tj = 7/T, tj = Q/T, for solid, dotted and dashed respectively 

sensitive prove of the viscosity. Unfortunately such constraints are not available at the 
moment. 

IX. CONCLUSIONS 

In this paper we have extended our studies in Ref. of the linearized hydrodynamic 
equations which describe the interaction of jets with quark gluon plasma. We have clarified 
the origin of the two hydrodynamic modes that we called sound and diffuson in Ref. ^ and 
related them to vortex free and rotational solutions. 

We then used the relativistic analogue potential flow to find the governing equations of 
these modes in an expanding background. For a static background, we subsequently used 
these equations to systematically study the interaction of a jet with the medium. Generally 
the jet can excite both the sound and the diffusion modes. The entropy produced by the 
jet-medium interaction fixes the strength of the diffusion mode relative to the sound mode. 
This entropy is localized in a narrow wake along the direction of the jet. 

The equations for the sound wave are given by Eqs. ()4.2|) . and ()5.3p . The unknown 
function in this formula, dF/dx, is related to the momentum transfered to the sound wave - 
see Eq. fl6.5|) . Similarly the flow fields in the wake are given by Eq. ()5.4p . and the amplitude 
A in this formula is fixed by the momentum transfered to the wake - see Eq. ()5.13|) . This 
momentum transfer is in turn related by Eq. (j5.14j) to the total entropy produced by the 
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jet. In summary, by specifying the total rate of energy loss and entropy production the flow 
fields at large distances are determined. 

Regardless of the mechanism of excitation, a supersonic jet in a static fluid leads to 
flow at the Mach angle. However, the different mechanisms of excitation may or may not 
lead to visible peaks in the final particle spectra. If the jet-medium interaction produces 
significant entropy, the wake contribution obscures the Mach angle in the final spectrum. 
Only very strong jet-medium interaction which does not produce significant entropy can 
yield hydrodynamic fields that give peaks in the azimuthal correlation functions at the Mach 
angle, A(/) ^ tt — 1.2 rad. For example, for an isentropic jet-medium interaction with a large 
energy loss dE/dx ~ 12 GeV/fm in a static medium, the resulting flow produces azimuthal 
distributions of the same order of magnitude as those observed by PHENIX and STAR 
0, 0|. Unfortunately, for moderate energy loss, dE/dx = 2 GeV/fm, the flow is unable 
to quahtatively reproduce the data. However, we argued previously that the expansion may 
amplify the sound wave and ultimately reduce the required dE/dx ji^. We have also shown 
that for the large dE/dx, the pt dependence of the correlation function has some of the 
qualitative features of the experimental measurements. 

If the jet-medium interaction does not produce significant entropy, the observation of 
peaks in the angular distribution is possible, but it is very sensitive to the characteristic 
scales involved. Systematic study of the viscosity, the jet-medium interaction time, the 
freeze out time, and the source size a, reveals that the prediction of a Mach cone from 
linearized hydrodynamics is quite fragile and depends on many of the microscopic details 
of the interaction of the jet with the medium. Further work is needed to clarify these 
microscopic details before a strong conclusion can be reached. 
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APPENDIX A: ABOUT THE EQUATIONS OF THE POTENTIAL 



1. Derivation of the dynamic equation 

In this section we derive a nonlinearequation for the potential introduced in section 4. 
As in the case of irrotational flow for ideal hydrodynamics, the only non trivial equation 
is the entropy continuity equation Eq. ()4.5j) the equation for the potential arises from re 
expressing Eq. ()4.5|) in terms of (j). Let us start by writing Eq. ()4.5|) as 

d,^d^<P+^d,d^<P = Q. (Al) 

As for a baryon free fluid the speed of sound can be expressed as 

we can express the first derivative on Eq. (jAljl in terms of derivatives of the temperature as 

Now we note that from the definition of (j) we obtain 

= d.d^cj) . (A4) 

Dropping the factor s/T and multiplying by T^, we obtain the following nonlinear equation 
equation for (f) 

4 - 1) d^ct>d,(t>d^'d''<p + d^(t>d^ct)d,d''<p = o . (as) 

For a baryon free fluid, the hydrodynamic equation reduces to the energy momentum con- 
servation. Thus, in the case where the initial conditions for ideal hydrodynamics do not 
introduce any vorticity, the equation derived is totally equivalent to solving the standard 
system of equation of conservation of the stress energy tensor. 



2. Bjorken solution as potential flow. 

As an example of potential flow let us look for a boost invariant solution of hydrodynamics 
with no transverse expansion [sH. We will assume that the initial conditions do not introduce 
any vorticity, and thus, all the properties of the flow field can be derived form a potential 
(f). As is a scalar, the requirement of boost invariance imposes that is a function of 
the proper time r = \Jt'^ — only. For this potential, the nonlinear differential equation 
Eq. ()A5|) takes the form 

( d V (\ , \ d \ 
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Using the relation between the temperature and the potential Eq. ()A4|) we obtain T as a 
function of r 

r = ^^M. (A7) 



Thus, Eq. flA6|l becomes a first order equation for the temperature 



-It + = , (A8) 

dr T 



which, by means of the equation of state leads to 



I d 1 , , , 



the solution of which is s oc 1/r, the well know Bjorken solutions [3l|. The velocity field 
can also be readily calculated from the potential From its definition, Tu^ = d^cf), we obtain 
that the the velocity field is = 1 or 

= l(t,0,0,z) , (AlO) 

T 



that coincides again with the result due to Bjorken [31[. 



APPENDIX B: FAR FIELD FROM SUMMING INFINITESIMAL DISTUR- 
BANCES 

We now want to connect our previous solution on ji| with our general expression Eq. ()5.3p 
for the perturbed fields in the ideal case. The goal is to show that both ways of calculating 
give the same expression in the region far from the fluid. However, extra terms from those 
in Eq. (j5.3j) appear that are only significant in the region close to the jet. 

As argued along the text, the potential flow describes sound waves, that are also described 
for the system of equations Eq. (j3.4p . Setting = 0, the system of equations ca be written 
as the wave equation for e. We will consider the solution for an infinitesimal disturbance 
that happens in an time interval dto around to that leads spherically symmetric disturbance 

edio(^o,x) = eo(Rfo) , (Bl) 
gdtoito,:)^) = goiRto) , (B2) 

where Uto = ^/{x + to)^ + for a jet moving at the speed of light along the —x direc- 
tion. The total disturbance is obtained by the addition of all the infinitesimal disturbances. 
The general solution for the spherically symmetric wave takes the form 

eto = ^ (/- (K - Csit - to)) + U (Rto + Csit - to))) (B3) 

The values of the two functions are set by the initial condition Eq. ()B1|) at to. As the 
incoming wave /_)_ is only important in the region of the order of the source size we can 
neglect it in the far field. 
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From the definition of tlie potential in tlie static case, tlie perturbed temperature can be 
as T' = dtf and tlius, by means of the equation of state, we obtain the potential 

^cito = ^{F{Rto-Csit-to))) , (B4) 

where F = c^f^/s where F denotes time derivative. The addition of all the infinitesimal 
disturbances leads to 

^= / dtoe{t-to) . (B5) 



oo 



RtO 



The solution Eq. ()B5|) can be rewritten as Eq. ()5.3p by introducing a 5— function as 

F {RtO -clit- to)) = r d^^F (-^^) 6 {RtO - c,{t - ^/v - to)) (B6) 

The integration of over to can be done now by solving the delta function. This is standard 
and after this solution, we can express 

5 {RtO - Csjt - ^/v - to)) ^ 6 {to - tj) 

Rdm c,((x + vt-0'-/5V)'^' 

with P defined in Eq. ()5.2|) . and t, is the solution of the 5— function. 

t,-t+^ = (-{^:, + yt-i)±^J{x + vt- - P^p^\ . (B8) 



V — c% 



The existence of solutions for the 5— function imposes a constraint in the integration such 
that the radicand in Eq. ()B8|) is greater than zero. The 6'— function in Eq. ()B5|) constraints 
further this integration. Imposing t > tj leads to two different integration regions: 

1. X + vt > c^g/v'^{x + vt — Pp), that includes the region inside of the Mach cone. In this 
case requiring that t > to leads to the field 

^ ( l'^+'"^~f^P p-vRt/c \ p 



2. X + vt < cl/v'^{x + vt — Pp), that corresponds to the field out of the cone The field is 
then 

ip = - ^ (BIO) 
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